An invariant distribution in static granular media 
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We have discovered an invariant distribution for local packing configurations in static granular 
media. This distribution holds in experiments for packing fractions covering most of the range 
from random loose packed to random close packed, for beads packed both in air and in water. 
Assuming only that there exist elementary cells in which the system volume is subdivided, we derive 
from statistical mechanics a distribution that is in accord with the observations. This universal 
distribution function for granular media is analogous to the Maxwell-Boltzmann distribution for 
molecular gasses. 
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PACS numbers: 45.70.-n Granular Systems 45.70.Cc Static sandpiles; Granular Compaction 81.05.Rm- 
Porous materials; granular materials 



Granular materials are complex systems character- 
ized by unusual static and dynamic properties. These 
systems are comprised of large numbers of dissipative 
macroscopic particles assembled into disordered struc- 
tures. The microscopic description of the system-state 
requires a very large number of variables. However, there 
is a very large number of microscopic configurations cor- 
responding to the same macroscopic properties. Edwards 
and coauthors p], Q have proposed that the complex- 
ity of static granular systems could be disentangled by 
means of a statistical mechanics approach reducing the 
description of the system state to a few parameters only 
@, i, B i, 0, H IE 03, HO, El- An essential part of 
Edwards' idea is that in static granular media volume 
plays the role held by energy in usual thermodynamics. 
Therefore, an understanding of the volume distribution 
function is the key to connect microscopic details of the 
system with macroscopic state variables. 

Since granular materials are dissipative, they can 
change their static configurations only when energy is 
injected into the system. The general idea underlying 
a statistical mechanics description is that the properties 
of the system do not depend on the kind of energy in- 
jections, but only on the portion of the configurational 
space that the system explores under such action. For in- 
stance, some classical experiments [13, |lj, |l5[ obtained 
different average packing fractions by tapping the con- 
tainer with different intensities and different numbers of 
times. Similarly, in an experiment by Schroter et al. |l6j |. 
reproducible average packing fractions were obtained by 
driving the system with periodic trains of flow pulses in a 
fluidized bed. More generally, various controlled pertur- 
bations (tapping, rotating, pouring, etc.) can produce 
packings with characteristic average packing fractions. 
Within a statistical mechanics framework, if the system 
volume (Vt) is the relevant state variable, the system 
properties should depend only on the achieved packing 
fraction and not on the preparation history. This hy- 
pothesis is examined in this paper using the largest sets 



of experimental data on particle positions presently avail- 
able. 

Experiments.- We analyze the structural properties of 
static granular packings produced in 18 different exper- 
iments, 6 with acrylic spheres in air and 12 with glass 
beads in water. The packing fractions p range from 0.56 
to 0.64. Three-dimensional density maps have been ob- 
tained for these systems using X-ray Computed Tomog- 
raphy [17]. Coordinates of the bead centers have been 
calculated for more than two million spheres with a pre- 
cision better than 1% of their diameters, which is better 
than the uncertainty arising from polydispersity (5% for 
the glass particles and 2% for the acrylic particles). 

The 12 samples of glass beads were prepared using the 
fluidized bed technique described in [Tf|; each sample 
consisted of about 145,000 spherical grains with diame- 
ters 250 ± 13^m placed in a cylindrical glass container 
with inner diameter 12.7 mm. The beads were fluidized 
with pulses of deionized water with the flow rate pulsed 
by a computer-controlled syringe pump. During each 
flow pulse the bed expanded until its height reached a 
stable value. After each flow pulse the bed was allowed 
to settle into a mechanically stable configuration. Pack- 
ing fractions in the range between 0.56 and 0.60 were 
obtained by tuning the water flux. The particle configu- 
rations have been studied for sub-sets of grains that were 
at a distance larger than four sphere diameters from the 
sample boundaries; there were about 90,000 grains per 
sample. 

The 6 samples with acrylic spheres were described 
in |l8[. In particular, two samples contained 140,000 
particles (90,000 in the internal volume) with diameter 
1.59 ± 0.02 mm in a cylindrical container with an inner 
diameter of 55 mm, filled to a height of 75 mm. The 
other four samples each had about 35,000 (15,000 in the 
internal volume) acrylic spheres with diameter 1.00±0.02 
mm in the same cylindrical container. 

Theory.- Since the original work by Edwards and co- 
authors, several statistical mechanics approaches have 
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FIG. 1: (a) Distributions of the Dclaunay cell volumes from 
18 experiments collapse onto a universal curve when plotted vs. 
(V — V m i„)/({V) — V m i n ). The inset shows the same distribu- 
tions plotted vs. V/cP. (b) The inverse cumulate distributions, 
F(V) = 1 — Jy f(V)dV), also show evidence of a collapse onto 
a universal curve when plotted vs. (V — V m i„) / ((V) — V m i„). The 
line is the inverse cumulate distributions for f(V,k = 1) (EqUJ. 
The inset shows the cumulate distributions plotted versus V/d 3 . 
The distributions are obtained from a statistical analysis over more 
than six million cells in 18 different experiments with three differ- 
ent kinds of beads in air and in water. The open symbols cor- 
respond to the 6 experiments with with dry acrylic beads [l8|l : 
circles sample A (p = 0.586); squares sample B (p = 0.596); stars 
sample C (p = 0.619); diamonds sample D (p = 0.626); triangles 
up sample E (p = 0.630); triangles down sample F (p = 0.64). 
The '+' corresponds to the 12 experiments with packing fractions 
0.56 < p < 0.60 made with glass spheres in of fluidized bed. 



been proposed to describe granular materials [l|, 0, 0, 
H, H, 0, H, EH • We make here the minimal number of 
assumptions in considering a granular system at rest that 
occupies a given volume Vt- We assume that the system 
can be subdivided into C elementary cells with volumes 
Vi (i — 1...C) with Xa=i v i = Vt- Under the sole as- 
sumption that any assembly of such elementary volumes 
will produce proper mechanically stable packings, we ob- 
tain the probability distribution for the cell-elementary 
volumes: 



X = (v) 



(2) 



where (v) — Vr/C is the average volume per elemen- 
tary cell, v m i n is the minimum volume attainable by the 
elementary cells and x 1S an intensive thermodynamic pa- 
rameter accounting for the exchange of volume between 
the elementary cell and the surrounding volume "reser- 
voir" . There is no need to specify the nature and kind of 
such "elementary cells"; the only assumptions are that 
they exist, that they can have any volume larger than 
Vmin, and that their assembly fills a volume Vt- Space 
can be divided arbitrarily into pieces. Common exam- 
ples of such partitions into space-filling blocks are the 
Delaunay and the Voronoi decompositions [2(J HH 
Such cells do not coincide with the elementary ones, but 
they might be assemblies of such elementary cells. For 
instance, in a Delaunay partition of a three-dimensional 
packing there are about six times more cells than in the 
Voronoi decomposition. The present theory applies to 
any degree of space partition made by any agglomerate 
local structure made of a given number k of elementary 
cells. From Eq.Q]the probability distribution for the vol- 
umes of such agglomerate of k elementary cells is 



f(V,k) 



(fc-1)! {{V)-V ml n) k 



exp 



V - V„ , 



(V) - Vmin 

(3) 

with (V) = Vt/C and V m in = kv m in being respectively 
the average and the minimum volumes for a given pack- 
ing. It follows from Eqs. [T] and [3] that 



X = 



(V) ~ Vm 



(4) 



p(v) 



-(v-v min )/x 



X 



(1) 



which is the average free-volume per elementary cell. 
Therefore, the intensive variable x is a measure of the 
kind and the degree of space-partition into elementary 
cells. 

We have obtained f(V, k) assuming that the cells are 
uniquely characterized by their volumes, and that any 
combination of C cells with arbitrary volumes Vi > v m m 
will produce a structurally stable, space-filling system of 
cells. This is possible in one-dimension only where, in- 
deed, the packing is an arbitrary disposition of grain cen- 
ters at distances larger or equal than v m in- In this case, 
the elementary cells are the Delaunay cells [22[ , which are 
the segments in between two successive grain centers; the 
distribution of their sizes is exactly described by p(v). In 
three dimensions, the Delaunay cells are tetrahedra with 
vertices on the centers of neighborhood grains chosen in 
a way that no other grains in the packing have centers 
within the circumsphere of each Delaunay tetrahedron. 
Clearly, in this case, the Delaunay cells are not uniquely 
described by their volumes and an arbitrary collection of 
cells is neither space-filling nor mechanically stable. 

Results.- We observe that the distributions of the De- 
launay volumes obtained from the 18 sets of data for 
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FIG. 2: The distributions of the Vorono'i cell volumes collapse 
onto a universal curve when plotted vs. (V — V m i n )/({V) — V m i n ). 
The theoretical line is f(V,k = 12) (Eq. [TJ- The inset shows the 
distributions plotted versus V/d? . The data and the symbols are 
the same as in Fig. [T] 



different kinds of beads in different media and different 
conditions are all described well by the same function, as 
shown Fig[lja). The collapse of the data onto a single 
curve was obtained using V m in = v / 2/12d 3 (with d the 
sphere diameter), which is the volume of a regular Delau- 
nay tetrahedron for four spheres in contact (the smallest 
compact tetrahedron [HI). The tail of the observed dis- 
tribution function for the Dclauncy volumes is exponen- 
tial, as in Eq. [Tjwith k ~ 1. The data fit well to the in- 
verse cumulate distribution, F(V) — 1 — J v . f(v', l)dv' , 
as FigOJb) shows; the agreement extends over four orders 
of magnitude and uses no adjustable parameters. The 
agreement suggests that Delaunay cells might be consid- 
ered candidates for the "elementary cells"; however, it 
was noted in [2l[ that, at small V — Vmin in the region 
where the spheres make contact (V < \/3d 3 /12), the em- 
pirical distributions deviate from the simple exponential 
form predicted by p(v), which indicates that other con- 
straints such as mechanical equilibrium should be also 
taken into account to correctly describe such region. 

An alternative way for dividing space into space-filling 
cells is the Vorono'i partition. In one dimension, the 
Vorono'i cell around grain 'V is constructed by taking 
the segment between the two mid-points between grains 
i — 1 and i, and between grains i and i + 1. The size 
of such a segment is equal to (u$ -I- Vi+±)/2 where Vi 
and Vi+i are the distances between the two consecutive 
couples of points. Therefore, the probability of find- 
ing a one-dimensional Vorono'i cell of size V is associ- 
ated with the probability of finding two successive De- 
launay cells with sizes v\ + «2 = 2V. This probability 

is fL'u™" 1 p( Vi )p( v ~ v i)dv\ oc (V - 2v m i n )exp[-(V - 
2v min J/\(V) - 2v mm )] ex f(V,k = 2), with (V) the av- 
erage size of the Vorono'i cell. In three dimensions, the 
Vorono'i cell can be also seen as the combination of sev- 
eral elementary cells with distribution f(V,k); however, 
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FIG. 3: The standard deviations of the Vorono'i volume distribu- 
tions in the 18 experiments fit the linear relation given by Eq. \E\ 
with k = 12. The dash line above (below) the solid line corresponds 
to k = 11 (k = 13). 

in this case the number of sub-cells involved is not fixed 
at k = 2 but depends on the kind of packing. 

We observe that data for over a million of Vorono'i cells 
from all 18 experiments for dry and wet packings of glass 
and acrylic spheres collapse to the the same distribution 
function (see Fig. ^ . In this collapse of the data there 
no adjustable paramaters, just (V) = Vr/ (number of 

grains) and V mm = b^/^/^J 2{2Q + 13\/5)d 3 ~ 0.694d 3 , 
which is the smallest Vorono'i cell that can be built in 
a equal-spheres packing [22| . Figure [2] shows that such 
universal distribution function is well described by Eq. [T] 
with k — 12, which indicates that about 12 elementary 
cells contribute in building each Vorono'i cell. Such a 
number is meaningful, since about 12 spheres are ex- 
pected to be found in the close neighboring of any given 
sphere in the packing. 

A further demonstration that such statistical distribu- 
tions are independent of the details of a sample and the 
method of sample preparation is provided by the behav- 
ior of the volume fluctuations, which can be calculated 
directly from the relation (see Eq. Q}, 

v 2 (V) = ((V-(V)f)= X 2 ?SYl , ( 5) 
which becomes, using d (V) jd\ = k (from Eq. H}, 

a(V) = Vk X = {V) ~J™ m . (6) 

The good correspondence of the data from the Vorono'i 
volume distributions with Eq. [6] (Fig. [3]) provides further 
evidence that (V) and k are the relevant control param- 
eters, independent of the system and its preparation. 

Conclusion.- We have shown that for a broad set 
of granular packings of different monodisperse spheri- 
cal grains, prepared by using different procedures, the 
local volumes are described by a universal distribution 
function (Eq. [3]) , which can be predicted from a minimal 
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set of assumptions. The only tunable parameter in the 
present theory is k, which has the value 12 for Voronoi 
decompositions throughout the accessible density ranges 
of the different static granular packings studied. Interest- 
ingly, for granular gasses the same empirical distribution 
(Eq. [3|) applies but with k = 5.586 [2H. [23|. A similar 
kind of distribution (Gamma distribution) with differ- 
ent values of the parameter k have also been observed 
in two dimensional Voronoi tessellation generated from 
disk packings [25j]. Therefore, the parameter k could be 
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